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The Renyi entropies of a massless Dirac fermion on a circle with chemical potential 
are calculated analytically at nonzero temperature by using the bosonization method. The 
O bosonization of a massive Dirac fermion to the sine-Gordon model lets us obtain the small 

mass corrections to the entropies. We numerically compute the Renyi entropies by putting 
> a massive fermion on the lattice and find agreement between the analytic and numerical 

^ results. In the presence of a mass gap, we show that corrections to Renyi and entanglement 

entropies in the limit m gap 3> T scale as e~ mgap / T . We also show that when there is ground 
state degeneracy in the gapless case, the limits m gap — > and T — > do not commute. 
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1 Introduction 

Entanglement entropy is a unifying theme in many different areas of theoretical physics 
today. In relativistic field theories, certain special kinds of entanglement entropy show 
monotonicity properties under renormalization group flow [TJ[2] . For conformal field theo- 
ries in (1 + l)-dimensions, numerical computation of the entanglement entropy provides a 
rapid way to calculate the central charge c. In the context of condensed matter physics, 
entanglement entropy can detect exotic phase transitions for systems lacking a local order 
parameter. The Ryu-Takayanagi proposal [3j|4] for computing the entanglement entropy 
holographically connects this circle of ideas to general relativity and string theory via the 
AdS/CFT correspondence [5}|7]. See refs. [8 12 for reviews. 



Recall the entanglement entropy is defined from a reduced density matrix p&. We start by 
partitioning the Hilbert space into pieces A and complement A = B. Typically A corresponds 
to a spatial region. We form the reduced density matrix pa = tr#p by tracing over the 



1 



degrees of freedom in B. Finally, the entanglement entropy is defined to be 



S = - tr p A log p A 



Related quantities are the Renyi entropies 



Sr, 



n 



logtr(p A ) 



(2) 



Note that the entanglement entropy can be determined from the limit S = lim n ^! S n . 

Given entanglement entropy's prominent role, it is surprisingly difficult to compute, even 
for free field theories in (1 + l)-dimensions. In this paper, we present some new results for 
the free, massive, Dirac fermion in (1 + l)-dimensions. We are particularly interested in 
thermal and finite size corrections to the entanglement entropy, and so we place our massive 
fermion on a torusj^] We allow for a nonzero chemical potential as well. 

Before proceeding further, let us briefly review the known results for the free fermion. 
Consider a massless Dirac fermion on the real line where the region A consists of p intervals 
whose endpoints are described by the pairs of numbers (x a ,y a ), a = 1, . . . ,p. In this case, 
the entanglement entropy takes the form (Tip 



S 



log 



Ua,b( X a - Vb) 



eP Ua<b( x a-x b )(y a -yb) 



c 



(3) 



where e is a UV cutoff and Co is a cutoff dependent constant. Given that the fermion is 
massless, we can use conformal symmetry to map the plane to a cylinder with either time or 



space compactified 15 . In the first case, we make the replacement (x — y) — > sinh irT(x — y) 



where T is the temperature. In the second case, we instead send (x — y) — > sin7r(:r — y)/L 
where L is the circumference of the spatial circle. 

If we turn on a mass m > 0, we lose conformal symmetry, and the computations get 
correspondingly more difficult. For the fermion on the real line at zero temperature and 
a single interval, the entanglement entropy can be expressed in terms of a solution to the 



Painleve V equation 14 . While in general only a numerical solution to the differential 



equation is available, one can find small and large mass expansions. For small mass, the 



^ prequel |13| to this paper considered the free massive scalar on a torus; the two papers can be read 
independently. 



2 See also ref. 15 
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leading log result is 

1, I 1 



5= - log - - -(mlXogml) 2 + 0( (ml) 2 \og ml) . (4) 
3 e 6 



At large mass, there is instead exponential suppression^] 



S^-J—e- 2m£ . (5) 
8 V ml K ' 

For multiple intervals, ref. [18] provides a small mass expansion. The leading log correction 
is instead — ^(mltlogme) 2 where It is the total length of all of the intervals in A. 

This paper contains two principal results. The first is a computation of the Renyi en- 
tropies for a region A consisting of multiple intervals where both 1/T and L are kept finite, 
the chemical potential \x can be different from zero, but m — 0. Previously, only the single 



interval Renyi entropy was available 19 20 . We compute the entanglement entropy from 
the Renyi entropies by analytic continuation. While the Renyi entropies are expressed com- 
pactly in terms of elliptic theta functions, our entanglement entropy is given as an infinite 
sum. In the limit where L — > oo or T — > 0, only the first few terms in the sum contribute, 
and we recover the cylinder version of eq. (|3j) . However, there is a subtlety in the T — > 
limit that we will return to shortly. 

The second principal result in this paper is a computation of the leading small mass 
correction to the entanglement entropy when both 1/T and L are finite. The correction is 



obtained using the equivalence of the massive Dirac fermion to the sine-Gordon model 21 22 
The result is expressed as a double integral over a product of elliptic theta functions. We are 
able to perform the double integral numerically and match the result to a numerical lattice 
computation of the entanglement entropy. We can also perform the integral in the limit 
1/T, L oo where we recover the multi- interval version of the leading log correction Q. 

We have found several interesting features of the entanglement entropy in the small 
mass and low temperature regime. The first is that the limits m — > and T — > of the 
entanglement entropy do not commute when there is ground state degeneracy. In particular, 
when the massless Dirac fermions have periodic boundary conditions around the circle, the 
ground state of the system is four-fold degenerate. If we first set m = and then take 
T — > 0, the system is not in a pure state. However, if we instead send T — > and then 
take m — > 0, the system will be pure. For pure states, the entanglement entropy of a region 
and its complement must be the same, S(A) = S(A). For example, consider the single 



s See also refs. 16 17 
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interval case, £ = x — y, of eq. (J3J) where we have conformally mapped to the cylinder 
(x — y) — > sin7r(x — y)/L. Then the entanglement entropy takes the form 

5=ilog(^sin^+c , (6) 

which clearly satisfies S(£) = S(L — £). If instead we take m — > first, we find a correction 
to eq. that reflects the ground state degeneracy. There is a similar correction in the 
multi-interval case. 

For Dirac fermions with antiperiodic boundary conditions, there is a unique ground state. 
The mass gap to the first excited state is m gap — tt/L. In this case, we can examine 
the corrections to the entanglement entropy in the low temperature limit T <C tt/L. The 
prequel |13| to this paper, based on an investigation of the massive (l+l)-dimensional scalar, 
conjectured that such corrections should be exponentially suppressed. Indeed, we are able 
to confirm this conjecture for the antiperiodic fermions: 

S(T) - S{0) ~ e" mgap/T . (7) 



Moreover, our result for the Renyi entropies allows us to determine the coefficient in front 
of the exponential factor. For the periodic fermions with m > 0, we find similar behavior, 
but in this case, our results come from a numerical lattice computation. 

The Lagrangian density for a Dirac fermion in what would be mostly minus signature for 
the metric if we had more than one spatial direction is 

C DF = - . (8) 

Our conventions for the gamma matrices are that {7^,7^} = 2rj fiu = 2(H — ). We define 
i|r = \j/t^/0_ ^y e choose gamma matrices 7 = a 1 and 7 1 = —ia 2 . 



2 Bosonization and Conformal Field Theory 



Following 14 , we consider the entanglement entropy of a free Dirac fermion on a torus with 



multiple intervals (u a ,v a ) (a — 1, ... ,p). Instead of having a single field on the n-covering 
space, we introduce n decoupled fields \P fe (k = — t ^,--- , rj ^) living on a single torus. 
They are multivalued around the branch points u a , v a around which they get phases e 1 ^ 
and e - * 2 ?, respectively. Here we define the single- valued field \l/ fc by introducing an external 



4 



gauge field fy k (x) = e^ x o dxl "- A ^ ( - x ^ k (x). It follows that the Lagrangian is given by 

C k = itf V(9 M + iA*)V k - m^ k . (9) 

The gauge field is almost pure gauge except at the branch points where delta function 
singularities are necessary to recover the correct phases of the multivalued fields 

OttIc P 

e^d v A\{x) = — [S {2) (x - u a ) - S®(x - v a )] . (10) 
0=1 

The partition function is then obtained in a factorized form 

n-l 
2 

Z[n}= ][ Z k , (11) 



h— — n ~' 
2 



where Z k is the partition function of the fc-th fermion coupled to the external gauge field A„ 

Z k = (jSAW**) , (12) 

with the current j£ = ^ k ^^ k . 

In (l + l)-dimensions, one can describe fermions in terms of non-local operators of scalars. 
The current is mapped to the derivative of a scalar field 

f k -> ^fMh , (13) 
and the Lagrangian of the k-th fermion becomes that of k-th real free massless scalar field 



0: C k = a-d^kd^^k- It follows from (12) and (13) that Z k can be written as the correlation 
function of the vertex operators 

p 

Z k = ([[V k (u a )V- k (v a )) , (14) 
0=1 

where the vertex operator V k is defined as V k (x) = e~ l n^ x \ 

The scalar field is a compactified boson with radius R = 2 so as to reproduce the partition 
function of a Dirac fermion on a torus]^] We have used the bosonization technique without 
specifying the spin structure of the fermion on a torus. We shall be more careful to distinguish 



4 In our conventions, R = v2 is the self-dual radius. 
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V 


sector 


{v x ,v 2 ) 


1 


(R,R) 


(0,0) 


2 


(R,NS) 


(o,§) 


3 


(NS,NS) 


(- -) 


4 


(NS,R) 


(1,0) 



Table 1: Conventions for fermion boundary conditions. 



the spin structures in the following. The torus is specified by two periods which we take 
as 1 and r = i/3 where (5 = 1/(TL) is the dimensionless inverse temperature]^] Let z be a 
holomorphic coordinate on the torus; then it has the periodicity z ~ z + 1 and z ~ z + r. 
The holomorphic part of the fermion on the torus satisfies four possible boundary conditions 



^{z + 1) = e 2mvi i)(z) , V(« + T ) = e 2lnU2 ^(z) 



(15) 



where i/j and is 2 take or |. The anti-holomorphic part satisfies the same boundary con- 
ditions as the holomorphic part. We denote the v = {y\,V2) sector where v = 1,2,3,4 
correspond to (0,0), (0, 1/2), (1/2, 1/2), (1/2,0), respectively (see Table [TJ. The correspond- 
ing partition function Z u is given by 



0„(O|t) 



7]{t) 



(16) 



Corresponding to the sector v, we can find a boson whose partition function agrees with Z u 
in the fermionic theory. The correlation function of the vertex operators on the torus in the 



is sector is given in ref. 23 



(O ei {zi, zi) • • ■ O eN (z N , z N )) u = Yl 



i<j 



d z A(0\r) 


2 




■&(zj - Zi\r) 




0„(O|t) 



(17) 



where R = 2 and O e is a vertex operator defined by O e (z, z) = e % ^ z,z ' . It follows from this 
formula that 



Zk,, 



Ii a <b^( u o. -u h \r) i?i(u -v b \r) 



Ua,bMUa -V b \T) 



(ed z M0\r)) p 



^(SE.(«.-*.)k) 



#„(0|r) 



3 We rescale the spacetime coordinates by L. 
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where we denote the partition function of the k-th fermion in the v sector by Z^ v . We 
normalize the partition function such that Z^ >v = 1 when there are no branch points. Since 
the theta function behaves as $i(z\t) ~ z in the small z limit, we put the UV cutoff e to split 



the coincident pointsjj Finally the total partition function (11) in the v sector is obtained 
as 



log Z„[n] 



n 2 -l 



log 



Ila<b$l( U a ~ U b\T) A(v a ~ V b \r) 



6n 

+ E 2 log 



n-l 
2 



M^E a (ua-v a )\r) 



^(0|r) 



(19) 



The Renyi entropy has the following form 

1 



1 — n 



(log Z>] -nlogZ.fi]) 



q i oM 
°n,0 ~r ° n l ■ 



(20) 



Here the first term is universal, 



S, 



71,0 



n+ 1 
6n 



log 



ria<b^l( U « - M fe! T ) - V b \ T ) 



Ua,bM u a ~ V b \T) 



(edM 



\r)) P 



and the second depends on the spin structure, 



n-l 
2 



^(jE.(«»-«.)|r) 



0„(O|r) 



(21) 



(22) 



Note that the Renyi entropies in the v — 1 sector are divergent because of the #i(0|r) inside 
the logarithm in S^\. We will have little to say about the massless u — 1 sector in what 
follows. 



2.1 Adding chemical potential 

In Lorentzian signature, a chemical potential is equivalent to introducing a constant time-like 
component of the vector potential A t = fi. In Euclidean signature, the chemical potential 
becomes pure imaginary, A tE = ifi. These considerations suggest that we can understand 

6 Here e is dimensionless. The dimensionful UV cutoff is eL. 
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the dependence of entanglement entropy on chemical potential by thinking about flat gauge 
connections on the torus. 

Let the vector potential be A k = at (He + ct x dx + . . . where a t and a x are constant 
and the ellipsis denotes terms responsible for the twisted boundary conditions around u a 
and t>&. Note that a t and a x are defined only up to gauge transformations which shift 
A k — > A k + 2im{\dtE + dx) where n is an integer. Such a flat connection contributes to the 



partition function through eq. (12). 



In the bosonized picture, the periodic scalar has boundary conditions along the thermal 
and spatial circles that are characterized by two winding numbers (w,w'): 



(z + l) = 4>(z) + 2tiRw 



[z + t) = (j)(z) + 2ttRw' . 



(23) 



The expectation value of the vertex operators is then computed by summing those over the 
topological sectors: 



e 3 z j) 



e (n°. 



j2i(/3a t w—a x w') 



(24) 



(w,w') 



When t = i(3, the (w, w') sector is related to the (0, 0) sector (see [23]): 



> z 3) 



(w,w') 



exp 



(0,0) 



2vri^, 







w' + Re(zj)w 



(25) 



From this result, we see that the effect of the flat gauge connection can be incorporated in 
the correlation function by making the shift 



e,- ej P 

3 3 



E 



(26) 



in the v dependent portion of the correlation function. 

Alternately, through the relation ^ k [x) = e^' x o dx '"' A ^^ x ^ k (x), we can trade the flat 
gauge connection for a shift in boundary conditions. From this expression, one may make 



the identifications a x = ±2tiv\ and a t = ±27ri/ 2 //3. Indeed, in eq. (17), we can rewrite the v 



S 



dependent term with the use of the formulae (88) in appendix |A| as 



It) 



0„(O|t) 



l9i( E i |^_^ r _ i/2 | r) 



7?l(-I/lT - U 2 \T) 



^1 



' Ei Gj^i 



+ 



2?r 



(27) 



To match the shift (26), we made the sign choices a x = 2-nui and a t = — 2iri , 2/f3. 



Introducing a chemical potential, by analytic continuation, is equivalent to introducing 
an imaginary a t = ifi. From the structure of eq. (26), it is clear that adding a chemical 



potential is also equivalent to introducing a real a x . This second equivalence makes it clear 
that the effect of chemical potential must be periodic with period 2n. Restoring dimensions, 
we see that the periodicity 2ir/L is precisely the energy level spacing on the torus. 



From eq. (27), one then obtains the partition function of the k-th fermion in the v sector 



with chemical potential 



k.v 



rL<b^iK» - u b \if3)§i(v a - v b \i(3) 



Ua,bM U a -V b \iP) 



(28) 



where £ t — L ^2 a {v a — u a ) is the total width of the intervals. The universal part of the Renyi 



entropy remains the same (21) and the part depending on the spin structure (22) is altered 
to 



°n,l 



n 



(29) 



Our result reduces to that of ref. 20 for one interval. 



2.2 Low temperature limit 

In this section, we consider a series expansion of Sn^ in the low temperature limit, r = i/3 — > 
zoo. We take advantage of the product representation of the theta functions (see appendix 



Bj) . The universal term S n fi on the right hand side of ( 20 ) becomes 



S. 



n,0 



n + 1 
6n 



log 



n a<b sin n(u a - u b ) sin ix(v a - v b ) 



lla,fe sin7r K -v b ) 



0{e 



-2-itIz 



(30) 
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in this limit. Note the entanglement entropy contribution can be straightforwardly recovered 
by setting n — 1, which in turn agrees with the spatial cylinder version of eq. (|3| reviewed 
in the introduction. 

However, to claim complete agreement with eq. (3), we need to check that does not 
contribute at zero temperature. Consider low temperature expansions of for the spin 
structures v = 2 and 3 corresponding to thermal boundary conditions. (The non-thermal 
spin structures v = 1 and 4 are given in appendix [Bj) For v = 2, defining r = ^2 a (v a — u a ), 
we find that 

s nl = Ss(n,r) + s 2 (n,r) , 



where 



n-l 
2 



5s(n, r) 



1 — n 



E 



cos 



n 



(31) 



(32) 



and 



s 2 (n,r 



1 — n 



E 



'— iy +1 1 / sm(irjr) 



j=antiperiodicl 

For v = 3, we find instead that 



j e 2 ^ -1 \sin(^) 



— 77, 



(33) 



c(3) 



2 ^ 



^— 1) J+1 / sin(7rjr) 



l-n^isinh7r/3j \sin(^) 



(34) 



Thus for spatially antiperiodic fermions, eq. (30) is the whole story at zero temperature 



while spatially periodic fermions get an extra correction 6s(n,r). 

To investigate the entanglement entropy, we take the n — > 1 limit. Much of this limit is 
straightforward: 



hm s 2 (n, r) = 2 > ; — — — e 



n— »1 



3=1 

oo 



hm = 2 £ 



j sinh7r/3j 
(— iy +1 1 — 7rjr cot(7rjr) 



3=1 



sinh 7r/?j 



(35) 
(36) 



These contributions vanish exponentially in the /3 — > oo limit. Our analysis of 5s(n,r) is 
incomplete. We find that 

5s(n,l) = 2 In 2 , (37) 



for all n, consistent with the fact that spatially periodic Dirac fermions have a ground state 
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degeneracy equal to four. For small r, we were able to obtain an asymptotic Euler-Maclaurin 
type expansion: 

(38) 



00 f9 2j ' — 1) 
lim 6s(n, r)=2j2 " ~ A — ~ B 2j C(2j) r 2j 



n— S>1 



J 



where Bj is a Bernoulli number and ((x) is the Riemann zeta function. Unfortunately, this 
expression is not Borel summable. 

2.3 High temperature expansion 

To investigate high temperature behavior, we use the modular transformation rules for the 
theta functions: 



$ x (z\t) = -{-ir)- xl2 e-^ z2lT '& 1 {z/T\-l/T) . 



(39) 



The modular transformations of the other theta functions are given in appendix |Xj The 
asymptotic form of the theta function depends on the value of z in the small limit: 



■& x {z/t\-1/t) = -2ie w sinh — + 0(e ? 



< z < 



(40) 



where r = i/3 was used. For 1/2 < z < 1, one may use the periodicity of the theta function: 



& x {z/t\-1/t) = e^ (2z - 1) ^ 1 ((l - z)/t\-1/t) 



(41) 



When Vp — Ui < 1/2, the leading term of the universal part S n> o of the Renyi entropy can be 
written 



5n,0 



1 + n) I irr z 



+ In 



n a<fe sinh J( "°- Mt) sinh 



P 



6n \ (3 



n a , fe sinh^pl 



(42) 



For v — 2 and 3, we find that 



D n,l 



'1 + n) irr 



6n /3 1 — n 



'-1)^ 1 /sinh^ 
J sinhf \s4nh^ 



n 



(43) 
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The entanglement entropy limit is given by 



lim S, 

n— yl 



11,1 



3/3" 



j'=i 



■1) 



^coth 



7rjr 



sinh^ 



(44) 



Similar results for v — 1 and 4 are given in appendix [Bj 

Note that the leading 7rr 2 //3 dependence cancels between S n $ and S^]. To recover the 
temporal cylinder version of eq. (|3]), we need to take (3 — > while keeping u a //3 and 
fixed. 



2.4 Mutual information 

The mutual Renyi information is an important measure of the entanglement between two 
intervals. Given two intervals A and B of length l\ and £2 separated by £3 on a circle of 
circumference L, the mutual Renyi information is 

I n (A, B) = S n (A) + 3 n {B) - S n (A U B) . (45) 

The definition makes clear that the mutual information is free of UV divergences, unlike the 



entanglement entropy. Using eq. (20), the mutual Renyi information of two intervals for a 



massless Dirac fermion on a circle at finite temperature becomes 

In(A,B) 



n 



1 



6n 



log 



M ei+ L e3 




r) 






r) 



71-1 

2 



n 



n-1 

2 



,9 (kh+l 2 
Uv \n L 


t)M0\t) 






r) 



(46) 



The logarithmic plots of the mutual Renyi informations for n = 2 in the u = 2,3,4 sectors 
are shown in Fig. [!} The mutual information is completely finite and positive. We let the 
width of two intervals A and B be i\ = £2 = L/10 and plot the mutual information with 
respect to the distance £3 between them. Since the two intervals are on a circle of radius L, 
I2 is symmetric under £ 3 — > L — £ x — £ 2 ~ £3 as is clear from the expression (46). The plots 
for n > 3 are qualitatively similar. 

We can use the high and low temperature expansions of Sn^ described above to get a 
better understanding of the behavior of I n . At large T, the theta functions can be replaced 



by hyperbolic sine functions, as in the expansion (42). In the v — 2 and 3 cases, expanding 



12 



0.2 



0.4 

h 
L 



0.6 



0.8 




-3 
-4 



0.0 



0.2 



0.4 

h 

L 



0.6 



Figure 1: The mutual Renyi informations of two intervals A and B of width t\ = ti = L/10 
with n = 2 in the z/ = 2 [Left], z/ = 3 [Middle] and v = 4 [Right] sectors. £3 is the distance 
between the two intervals. The blue dashed and orange solid curves are for /3 = 10, 1/5, 
respectively. 



the hyperbolic sines, for £ 3 < L/2 we find 7 2 ~ e 27rf 3^ w j 1 ji e f or £ 3 > L/2, by symmetry, 
I2 ~ e -27r(i-<i-4— ^s)r_ fj^g = 4 S pi n structure, however, develops an order one contribution 



to the entanglement entropy at high T, as can be seen from the expansion (97). At low T 



the theta functions are replaced by sine functions, as in the expansion (30). However, for 
v = 2, there is an extra contribution from 5s(n,r) because of the ground state degeneracy. 
Before closing this section, we compare our findings to the holographic computation 



24 , 25 where the mutual information undergoes a phase transition as the distance between 
the two intervals increases, i.e., I(A,B) 7^ for small £3 while I(A,B) = for large £3. In 
our case, the finite volume and finite number of degrees of freedom prevent a phase transition 
from happening. However, for large temperatures the mutual information exponentially falls 
off as £3 is increased for the "physical" (R,NS) and (NS,NS) fermions. 



3 Bosonization and the Sine-Gordon Model 

We used the bosonization technique to compute the entanglement entropy of a free massless 
Dirac fermion. Even after turning on the mass, one can still employ the bosonization from 
massive Dirac fermions to the sine-Gordon model: 



C SG = ^-0^8^ + X cos <f) , (47) 

S7T 



where A is proportional to the mass of the Dirac fermion: A = 26,27. Then, the 



leading correction of the partition function Z^ v starts from the 0(A 2 ) term due to the 
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charge conservation of vertex operators: 

Z Kv {m) = Z Kv + - / d 2 xd 2 y (cos <f>{x) cos <f>{y) Y[V k {u a )V. k {v a )) v + 0{X i ) . (48) 

J 0=1 



The integrand can be evaluated by using eq. (17) as follows 



p 1 p 

(cos0(x) cos(f)(y)Y[V k (u a )V- k (v a )) u = -(Y[V k (u a )\L k (v a ))u [A kt „(x,y) + A k)V (y,x)] 



a=l 



a=l 



(49) 



where 



A kjV (x,y) = 



• fe E<,K-«°) 



+ x-y\r) e<9 2 tfi(0|r) 



n 

a=l 





- x |t)#i(u 


- 


1?l(«a 


- x|r)^i(v a 


-?/l r ) 



(50) 



At leading order, the Renyi entropy is given by 

S£0(m) = 5^(0) + C„m 2 + 0(m 4 ) 
where the coefficient C n of m 2 is defined by 



(51) 



II — L 

i i r x 2 > 

Cn = " A 21 T\2 I d2xd2 V Y] A k;U (x,y) 

1 — n47r 2 (eL) 2 y ,. „_i 



(52) 



The four dimensional integral is too complicated to evaluate analytically and we shall rely on 
a numerical computation after isolating and showing the trivial nature of the UV divergence. 
Since 2k/n < 1, there are no poles at x,y — u a ,v a in the integrand A k)U (x,y). A possible 
divergence comes from the point x = y where $i(y — x\t) ~ y — x. Expanding the remainder 



of eq. (50) and summing it over k, one obtains the following series around x = y: 



n-l 
2 

E 



n-l 
2 



U kJ2A : a - Ua) + * - y\r) 



U kTia{ T Ua) \r) 



n 

a=l 





- a:|r)^i(« a 


-Z/I r ) 


#1 («o 


- s|r)i?i(v a 





(53) 



1 + 0((x - y) 2 , (x - y) 2 ,\x - y\ 2 ) 
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Therefore, the singular part of the integrand is 



n — l 
2 



^2 A ki „(x,y) 



\x - y\ 



[l + 0((x - y) 2 , (x - y) 2 , \x - y\ 



(54) 



The integration measure gives a factor of \x — y \ near x ~ y, and we end up with a single pole 
there. This single pole gives rise to the UV divergence after integration, but the divergence 
is independent of the size of the intervals. Since we are interested in the physics depending 
on the size, we will throw the divergence away and get a finite result in the end. Fig. [2] shows 
the result of a numerical integration of eq. (52) for one interval of width v± — u\ — t/L. We 



find good agreement with a lattice computation described in section |4| 




t 
L 



Figure 2: The t dependence of the 0(m 2 ) correction to the n = 2 Renyi entropy for v — 2. 



The curves are produced by numerical integration of (52). The points are from a lattice 
computation. From top to bottom, /3 = 1/2, 1, and 2. 



The IR divergence is absent on a torus, but it appears in the flat spacetime limit. It 
is worth looking into what happens in this case. The function Ak )V (x,y) reduces to the 
correlation function of the vertex operators on a flat spaceQ 



x 



n 



a=i 



(v a - x)(u a - y) 



(u a - x)(v a - y) 



(55) 



The most divergent term will come from the region x,y ~ A where A is the IR cut-off scale. 



7 The variables x, y, u a , Vb and e are dimensionless on a torus, but they have dimensions of length in the 
flat spacetime limit. We will use the same symbols for simplicity. 
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The expansion of A kjU (x,y) around large x and y is enough to compute the IR divergence: 



71-1 

2 



^2 A kA X 'V) 



n-1 
2 



\y - x\ 



n + 



n 2 — 1 [\x 


2 (y + y) - 




y 


2 (x + x)] 2 


Yin 


2 


X 


4 


y 


4 



(56) 



where i t = ^2 a {v a — u a ) is the total length of the intervals. The leading term is an £ t 
independent IR divergence, and we drop it below. Performing the integral over x and y, we 
obtain 



1 -I- n 

S n {m) = S n (0) - -^(rn£ t ) 2 log 2 A + 



(57) 



This small mass expansion is strikingly similar to the result ^ of refs. 14,18 reviewed in 
the introduction. 

We are working in a limit m 1/L,T, and our IR cutoff is naively given by the size of 
the torus L and /3. If we can commute the order of limits, we may identify the IR cutoff 



instead with the inverse mass of the Dirac fermion, A = 1/m, and then our result (57) agrees 



with Q. We will see below that the limits commute for the v = 3 spin structure but not for 



v — 2. In the v — 2 case, there is an extra contribution from 5s(n, r) in eq. (31) that needs 
to be removed when the limits are exchanged. 



4 Massive Fermion on the Lattice 



The Hamiltonian of a Dirac fermion on a circle of radius L can be derived from the corre- 
sponding Lagrangian density (J8l: 



(5f 



H=f <M t (-i7°7 1 3r + ra7°)* • 
Jo 



To put the fermion on the lattice, we discretize the circle into points with a lattice 
separation e = L/N 



N 

i=i 



,t ^3 



(59) 



16 



The canonical anti-commutation relations are {^j, a , ^\ g] — $jk 5 a p, and a, (5 are the spinor 
indices. To diagonalize the Hamiltonian, we expand the Dirac field as follows: 



N-V! 



To satisfy periodic {v\ = 0) or antiperiodic (y\ = 1/2) boundary conditions around the circle, 
we set 6 a = The energy u(9 a ) is defined by 



uj{6 a f = m 2 + ™^ • (61) 



This dispersion relation exhibits the classic doubling problem of fermions on the lattice. Our 
concern with finite size effects, however, introduces an additional subtlety. If we take N 
even, then we get two copies of either a V\ = or a v\ = 1/2 fermion. If we take N odd, then 
the second copy has the continuum spectrum with spatial periodicity opposite that indicated 
by v\. Numerically, we have observed that the entropy in this case corresponds to a v\ = 
plus a V\ = 1/2 fermion. To keep things simple, we will assume N is even from now on and 
then divide our entropies by two when comparing with the analytic results from earlier in 
the paper. 

The u a and v a are normalized such thalF] {b(9 a ), = S a b and {d(9 a ), S(9b)} = 6 a b- 



The vectors u a and v a satisfy the discretized Dirac equations 

, „ , n sin 9 n 1 \ 
uj{8 a ) 1 + —^ 1 1 -m\u a = 0, 



L i9a )^ + ^.^ +m \ Va = . 



(62) 



We demand that the u a and v a satisfy the normalization and orthogonality conditions 

u l u a = v\v a = 2u(6 a ) , 
u\v N - a = v\uN- a = • (63) 



One can explicitly find the vectors satisfying ( 62 1 and ( 63 1 



u„ 



/ , n . sin 8 a / , n . smi 



e 



t * . smf a / rn s suit 
v a = [ -\/u>(e a ) -i, ^uj{6 a ) + — 
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The Hamiltonian is diagonalized ag^] 



N-vi 



h= J2 uj(e a )[b\e a )b(e a ) + d\e a )d(e a )\ . (64) 

a=l— v\ 

In the lattice model, the fermion number operator is given by 

N-ui 

F = E {b\e a )b{e a )-d}{e a )d{e a )) . (65) 

a=l— i/i 

Introducing a chemical potential \x conjugate to F, the density matrix can be written in 
terms of H in the standard way: 

(_l\{l-^ 2 )F e -{H+nF)/T 

P = ti[(-l)( 1 - 2 ^) F e'( H +^ F )/ T ] ' ^ 

We have introduced a factor of (— 1) F to allow for spin structures periodic in the time 
direction. Expectation values are defined as (X) = tr(pX). A short calculation yields the 
two-point correlation function of two \& fields: 

to .*t v = J_ e «.y-») [ ( i + sinh (M ) 

\ 3 kl 2L J_,^ y cosh(/3/i) + (-1)2-2+1 co8k{pu)(6 a )) J 



sm6 a m \ sinh(/3a;(6> a )) 

+ ,n \ a 3 H 77TT Cr l 



W (0 a ) e J w(0 a ) V (-l) 2 ^+ 1 cosh(/3/i) + cosh(/3a;(^)) 



• (67) 



Note that the argument in section 2.1 implies the 1/2 = sector is obtained from the 1/2 = 1/2 
sector by shifting p — > p — in/ (3. The form of the two-point function is consistent with this 
observation. 

It is possible to calculate Renyi entropies from the matrix C^ v > = e(\l/\Ed). Consider a 
region A, which may consist of many disjoint subintervals of the circle, and the corresponding 



reduced density matrix p&. We restrict C-^ such that j and k run only over sites in A. Call 



the restricted two-point function C A U \ Remarkably, for a free spinor field, the reduced density 
matrix p^ ~ e~ HA can be written in terms of a free particle Hamiltonian Ha = £hb)jOk 
(see for example [9|[l0] ) . Moreover, there is a one-to-one correspondence between eigenvalues 



9 We remove the infinite constant coming from the commutation of d and <v . In other words, we fill out 
the Dirac sea. 
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Xj of Ca and the energies ef 



A, 



1 



1 - (-l) 2 ^ 

Given this relation, it is a short exercise to demonstrate that the Renyi entropies are 

1 



n 



tr log 



c 



(68) 



(69) 



where Cjk is the restricted two-point function. 



Before proceeding, we make two quick observations about the eigenvalue distribution of 



Cjf\ From the trace of C^' , we see that when the chemical potential vanishes, ^ - Xj 



where n is the total length of A. Next, from the relation (68), it is clear that in the thermal 
case (b> 2 = 1/2), Xj is bounded between zero and one. Provided < Xj < 1, we can take a 



sensible n — > 1 limit of (69) and derive the entanglement entropy 



S (u) 



tr 



c£°)]og(l + log 



(70) 



4.1 Comparison to the analytic results 

In section |2j using bosonization, we obtained analytic formulae for the Renyi entropies 
of a massless Dirac fermion. To gain confidence in our methods, we compare the lattice 
calculation for a massless fermion with these analytic formulae. Consider the Renyi entropy 
for n = 2 of two intervals of width t\ and ii separated by a distance of £3. In Fig.|3j we plot 
the entropies in the v = 2, 3, 4 sectors by changing the distance £3 with fixed l\ = £2 = L/10. 
The blue (dotted) and orange (solid) curves are the analytic results for (3 = 1/5 and 10, 
respectively. The dots are plotted using the lattice computation, which nicely agree with the 
analytic curves up to a constant. Since the Renyi entropy is always UV divergent, a constant 
shift is allowed to match the analytic and numerical results. The other case is studied in 
Fig.|4]by varying £2 with fixed £\ = £3 = L/10. The analytic and numerical results perfectly 
fit each other for various temperatures. 
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Figure 3: The Renyi entropy for n = 2 of two intervals of width £\ = £ 2 = L/10 whose 
distance is £ 3 . The v — 2 [Left], v = 3 [Middle] and v — 4 [Right] sectors are depicted. The 
curves are analytic and the dots are numerical. The blue dotted and orange solid curves are 
for /3 = 1/5 and 10, respectively. 




Figure 4: The Renyi entropy for n = 2 of two intervals of width l\ = L/10 and l<i- The 
distance between the intervals is fixed to £3 = L/10 and £2 is varied. The v = 2 [Left], v = 3 
[Middle] and v = 4 [Right] sectors are depicted. The curves are analytic and the dots are 
numerical. The orange solid, blue dotted and black dashed curves are for = 1/10, 1/5 and 
1. 
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4.2 Small mass vs. small temperature 



An interesting feature of the (R,NS) fermion is that the T — > and m — > limits do not 
commute. In the zero mass and zero chemical potential limit, the two point function becomes 



lim C jk 



N-vi 

/v e 



2iV 



0=1— fl 



l + a 3 sgn(sin a ) ^tanh 



4^2 — 1 



(71) 



If we further take the zero T limit of (71) for the (R,NS) and (NS,NS) fermions, we obtain 
respectively 



lim lim C 



(2) 



1 



T-+0 m-s>0 



'jk 



103 

N 



. 7r(j - fe) 
sin 



lim lim C 



(3) 



T-5-0 m-s>0 



^ jfc + (1 -Sjk) 1 ^ 



. 7r(j - fc) 
sin 



cot 



CSC 



Because the (NS,NS) theory is gapped even for m = 0, we find 



(72) 
(73) 



(3) 



lim lim C 4h 

T^0m->0 ■> 



(3) 



lim lim C„- t 

m->OT^0 J 



However, in the (R,NS) case, we find instead that 



(2) 



lim lim 

m->0T->0 3 



lim lim CfJ + ^7 



COS 



7r(j - fc) 



(74) 



Let us restrict to the case where A is a single interval of length n. It turns out that 
linv^o linij^o = (7r and lim^o linim->o Cjfe = Civs have the same eigenvalue spectrum, 
provided n is even. There is a similarity transform which relates the two 



C NS -M = M-C R 



(75) 



where 



M. 



r 0' 1)tT . (j - 1)7T 

cos — On+l-j,k Vi sin 



(76) 



AT " T1 ^" ' j\r 

^ < j,k < n. (For odd n, the eigenvalue spectra must then approach each other in the large 
N limit by continuity.) This equivalence means we can compute [lirnr->o> liiV-»o] ^(A) for 
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Figure 5: The single interval Renyi entropy for v = 2: n = 2 [Left] and n = 3 [Right]. In 
both cases, the \im T ^ lim m _;. curve (orange) lies above and the lim m ^ lim T _;.o curve (blue) 
lies below. The points were computed using the lattice. 

the (R,NS) fermion using our bosonization results: 



lim, lim 

_T->0 m->-0. 



S?(m,T) = W5f)(0,T) - Sf(0,T)) 

T— >0 



lim (S„ I S~ i 
Ss(n, r) . 



(3)> 



(77) 



For the entanglement entropy, we find that [limj^Oj lim m ^o] S^(m, T) = 5s(l,r). The non- 



commuting nature of these limits is shown in Fig. |5| Numerics suggest that the result (77) 
holds for multiple intervals as well. 



4.3 Small mass and temperature 



For theories with a mass gap m gap , ref. [13] conjectured that the temperature dependent 
portion of the entanglement entropy should have an exponential scaling dependence e -mgap / T 
in the range m gap ^> T. More precisely, given an interval A and its complement B, the 
conjecture posits that 



S A (T) - S B (T) ~ S A (T) - S A (0) ~ e~ m ^l T 



In this section, we provide further evidence for this conjecture. 

For the u\ = 1/2 fermions, the ground state is gapped with m gap = ir/L. Reassuringly, 
our low temperature expansions (36) and (94) for lim n ^i and lim n _;.i yield precisely 
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such scaling behavior, and we get the prefactor: 



S A (T) - S A (0) = ±4(1 - Trr cot(vrr))e-" /LT + 0(e- 2n/LT ) , 
S A {T) - S B (T) = T 47rcot(vrr)e-^ LT + 0( e - 27T/LT ) , 



(78) 
(79) 



where the top choice of sign corresponds to v — 3 and the bottom to v — 4. The region A 
is taken to have size rL. Similar scaling behavior holds for the Renyi entropies and can be 



computed from eqs. (34) and (92). 



We also investigate the scaling behavior for spatially periodic v% = fermions where we 
introduce an m ^ by hand. In this case, we have no analytic results to offer, but we can use 
the lattice to calculate the entanglement entropy numerically. We compute 5S = S(T) — S(0) 
for the (R,NS) fermion and a single interval. Fig. [6] clearly shows e~ m//T scaling in the region 
m^> T, both for small mass mL = 1/10 and large mass mL = 10. 



log(tfS) 




in 
T 



in 
T 



Figure 6: The entanglement entropy difference 5S = S(T) — 5(0) for (R,NS) fermions: [Left] 
mL = 1/10; [Right] mL = 10. The points are computed from a lattice, and the lines are fits 
with slope -1. From bottom to top, E/L = 1/10, 3/10, 1/2, 7/10. 



5 Discussion 

Our bosonization method of computing the Renyi entropy for a massive Dirac fermion is 
perturbative in the mass, and we would like to do better. As reviewed in the introduction, 



in flat spacetime, ref. 14 obtained a non-perturbative relation between the single interval 



Renyi entropy and a solution to the Painleve V equation. This non-perturbative relation 



uses a result of ref. 28 for the sine-Gordon model. Similar arguments may be useful for 
investigating the behavior of the massive fermion on a torus. 

Another possible non-perturbative approach is to use the lattice. The two fermion cor- 
relation function matrix C A that we derived above is Toeplitz. The Renyi entropy can be 
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expressed in terms of a contour integral over the characteristic polynomial of C A U \ Math- 
ematical techniques such as the Szego limit theorem and further generalizations such as 
the Fisher-Hartwig formula are available for taking such determinants. Indeed, such tech- 



niques have already been used to study the XY model 29 30 . Through a Jordan- Wigner 
transformation, the continuum limit of the XY model can be related to nonrelativistic free 
fermions. 

A field theory with a mass gap can be implemented geometrically by putting a gauge 
theory on a compact space. Such field theories sometimes have holographic duals with 



AdS geometries where the compact space is the conformal boundary 31 . Several authors 



have studied holographic entanglement entropies in these backgrounds 32 -34 . In the case 
of the mutual information, for strip like regions, there is a "phase transition" where the 
entanglement entropy is nonzero for two strips close together but vanishes once the strips 
become sufficiently far apart. In our case, we do not expect to have a phase transition given 
that we have a finite number of degrees of freedom and work in finite volume. Nevertheless, 
we do see that the mutual information is exponentially suppressed for large separations 
and high temperatures (at least for the "physical" v — 2 and 3 spin sectors). In the case of 
temperature dependence of the entanglement entropy, holographic examples typically predict 
that quantities such as Sa(T) — 5^(0) and Sa(T) — Sa{T) vanish exactly when T <C m ga _ p . In 
our case, we again see instead exponential suppression. Holographic theories are supposed to 
describe strongly coupled large- N field theories and the large- A r effect can drive the system 
to a phase transition. Presumably, we would need 1/N corrections to see holographically the 
exponential behavior observed in this paper. Perhaps these 1/N corrections could be studied 
by introducing higher derivative corrections, additional saddle-points in the path integral, 
or non-perturbative objects such as D-branes and orientifold planes. 

Another interesting direction for future study is to introduce interactions between the 



fermions. It is well known 21,22 that the sine-Gordon model, for more general choice of the 
interaction parameter A, fermionizes to the Thirring model which has a quartic interaction 
term. On the one hand, such a quartic interaction is not compatible with the replica trick 
where we replaced a single fermion field on the n-covering space with n decoupled fields 
living on a single torus. On the other, one could certainly use bosonization to treat the 
quartic interaction perturbatively. 
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A Theta Function Identities 

oo 

M^\r) = 2e 7riT/4 sin(7rz) -g m )(l -yq m ){l -y^q m ) , (80) 

m=l 

oo 

$ 2 (z\t) = 2e ntT/4 cos(nz)Y[(l-q m )(l + yq m )(l + y- 1 q m ) , (81) 

m=l 

oo 

^( Z \r)=l[(l-q m )(l + yq m - 1 / 2 )(l + y- 1 q m - 1 / 2 ) , (82) 

m=l 

oo 

Mz\t)= n( 1 -? m )(l-^ m " 1/2 )(l-Z/" 1 ? m_1/2 ) , (83) 

m=l 

where y = e 2mz and q = e 2mT . We also have the S-duality relations 





\r) 




{-ir)- 1/2 e- niz2/T A(z/r\ 


~ 1/r) , 


(84) 


d 2 {z 


\r) 


= (" 


-try^e-^^Mz/rl - 




(85) 




\r) 


= (" 


- ir )-V2 e -*i* 2 /T0 3 (3/ r | _ 


1/r), 


(86) 




\r) 


= (" 


-ir)- l / 2 e- mz2/T $ 2 (z/T\ - 


1/r). 


(87) 



The periodicities of the elliptic theta functions yield 



2 (z|t) = -#i(z-1/2|t) , (88) 
Uz\r) = -y-WqWM* ~ V 2 - r/2|r) , (89) 
Mz\r) = ty' 1/2 q 1/8 A(z ~ t/2\t) . (90) 
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B Time Periodic Spin Structures 

In the v — 1 sector, we regulate by introducing a small chemical potential fi <C 1/(3. 
The large (3 expansions for v — 1 and 4 are 



n— ^1 



2 



n,:L 1 - n ^ , 



7rkr 


+ 2 log 


up 


sin 


n 




2 


1 


( sin(7rj 


r) 



J'=l 



,(4) = 2 1 / sin(7 

171,1 1 - n j sinh vr/3j ^ sin ( 



sin(7rjr" 



— n 



n 



lim S^l = lim 

n-»-l ' n->l 1 — /j, 



2 



7rA;r 



sin 



+ 2 log 



fe^O 



oo 1 1 

4 7 e ^,- _ i ( 1 ~ ^ r cot(Trjr)) , 
j'=i 7 



- 2 E 



11 — 7rjr cot(7rjr) 



^ j sinh 7r/3j 
The small /3 expansions for v — 1 and 4 are 



(91) 
(92) 



(93) 
(94) 



c(l) _ 



lim S\ { 



1 + nnr 2 2 
~&n T + 



1 — n 



n-l 
2 



sinh 



7rA;r 



ra/3 



sinh ( ?f 



-je^-l ^ sinh (g) 



— n 



(1) _ tit' 



- 4 E 



1 1 



3/3 ^ j e^HP - 1 



2 



+ lim log 



sinh 
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+ 2 log 



irjr 
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+ 2 log 
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(95) 
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q(4) 



ra-»l ' 



1 + n nr 2 2 
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nkr 
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